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O ' Abstract 



Understanding the mechanism of electroweak symmetry breaking and the origin 

of boson and fermion masses is among the most pressing questions raised in contem- 

f^ i porary particle physics. If these issues involve one (several) Higgs boson(s), a precise 

1^ ■ measurement of all its (their) properties will be of prime importance. Among those, 

the Higgs coupling to matter fcrmions (the Yukawa coupling). At a Linear Collider, 
the process e~^e~ — > iiH will allow in principle a direct measurement of the top-Higgs 
Yukawa coupling. We present a realistic feasibility study of the measurement in the 
5_^ ■ context of the TESLA collider. Four channels are studied and the analysis is repeated 

C^ I for several Higgs mass values within the range 120 GeV/c^ - 200 GeV/c^. 
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1 Introduction 

The gauge sector of electroweak interactions has been checked to coincide with the Stan- 
dard Model (SM) prediction to the per-niil level, at LEP and SLC. On the contrary, there 
is no direct experimental evidence for the Higgs mechanism, supposed to be responsible 
for electroweak symmetry breaking and the generation of masses. Direct search of the 
Higgs boson at LEP yields the lower hmit [1]: Mh > 114.4 GeV/c^ at 95% CL. Precision 
measurements on the other hand give [2]: Mjj < 240 GeV/c^ at 99% CL. Once a Higgs 
particle is found, if ever, all its properties should be measured precisely to completely char- 
acterise the Higgs mechanism. Among those, the coupling of the Higgs boson to fermions 
(the Yukawa coupling), which is supposed to scale with the fermion mass: 

9ffH = ^ (1-1) 

where QffH is the Yukawa coupling of a fermion f of mass nif and v is the vacuum 
expectation value of the Higgs field, v = (\/2Gf)~^/^ ^ 246 GeV. 

The top quark is the heaviest fermion, thus the top-Higgs Yukawa coupling should be 
the easiest to measure. If Mh > 2 * mt, this parameter can be measured through the 
branching ratio of the Higgs boson decay into a pair of top quarks. Otherwise, i.e. for 
lower values of the Higgs boson mass, the process e'^e~ -^ tiH allows in principle a direct 
measurement of this coupling. 

Feasibility studies of the measurement of the top-Higgs Yukawa coupling via the process 
e^e~ -^ tiH at a Linear Collider have already been performed [3] [4] for a Higgs boson 
mass of 120-130 GeV/c^. This is the most favourable case (taking into account the lower 
mass bound) as the cross-section of this process decreases with increasing Higgs boson 
mass and as a Higgs boson of such a mass decays predominantly to a pair of b quarks, 
allowing a very effective signal and background separation using b-tagging algorithms. One 
of the studies ([4]) showed that a neural network analysis was essential to get a precise 
result. We repeated this work and extended it up to Mh = 150 GeV/c^. When Mh ^ 
135 GeV/c^, the H -^ W^W^ decay mode dominates. This channel was also studied, for 
masses up to 200 GeV/c^. 

2 The process e+e -^ tiH 

The lowest order Feynman diagrams contributing to the e~^e~ -^ tiH process are shown 
in figure 1. The amplitude of the diagram where the Higgs boson is radiated from the Z 
boson is not expressing the top-Higgs Yukawa coupling. However, since it modifies only 
slightly the cross-section of the process, it can safely be neglected. The cross-section and 
the top-Higgs Yukawa coupling thus verify to a good approximation: a^+^-^^jn oc gf^n- 

For this work, the following assumptions were made: mt = 175 GeV/c^ and BR{t -^ 
Wb) = 100%. The Higgs branching ratios were calculated with the HDECAY [5] program. 






Figure 1: Lowest order Feynman diagrams of the process e^e — > ttH . 

The values obtained for the H ^ bb and H -^ W~^W~ modes, which are the main decays 
within the Higgs mass range considered in this paper, are shown in table 1 and figure 2. 



Mh (GeV/c2) 


BR{H -^ bb) 


BR{H -^ W+W-) 


a(fb) 
(^ = 800 GeV) 


120 


67% 


13% 


2.50 


130 


51% 


30% 


2.17 


140 


33% 


50% 


1.88 


150 


16% 


70% 


1.64 


160 


3.1% 


92% 


1.44 


170 


0.76% 


97% 


1.25 


180 


0.48% 


93% 


1.09 


200 


0.23% 


73% 


0.80 



Table 1: Higgs branching ratios for the H ^ bb and H — > VF"*"!^" modes (as given by 
HDECAY) and cross- section at lowest order of the process e^e" — > iiH (as given by 
CompHEP), for various Higgs mass values and for y/s = 800 GeV. In the calculation of 
the cross- section, initial state radiation and beamstrahlung were taken into account. 



Previous studies showed that a center-of-mass energy close to the kinematical limit of 
the iiH production process does not allow to measure the coupling with a satisfactory 
precision, due to the tiny cross-section. Center-of-mass energies of 800 GeV and 1 TeV 
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Figure 2: Higgs branching ratios (dashed lines) for the H ^ bb and H -^ W~^W~ modes 
(as given by HDECAY) and cross-section (solid line) at lowest order of the process e~^e~ -^ 
tiH (as given by CompHEP), for various Higgs mass values and for y^ = 800 Ge V. In 
the calculation of the cross- section, initial state radiation and beamstrahlung were taken 
into account. 



have shown to be far better for the measurement. A run at 800 GeV being a very hkely 
possibihty for the Linear Colhder, this wih be the default value of our study. The lowest 
order cross-section of the process (as given by the CompHEP [6] software) is shown in 
table 1 and figure 2 for the Higgs mass range considered and a center-of-mass energy of 
800 GeV. 

The 0{as) corrections to the e^e~ -^ ttH process have been calculated by several 
groups [7]. At a center-of-mass energy of 800 GeV, they affect the total cross-section by 
less than 5% and were thus neglected in this study. 



3 Measurement of gun 

For a particular analysis yielding a selection efficiency of the signal e^gf "" and a purity of 



the selected sample 



sample 



and assuming an integrated luminosity L, the statistical and 



systematic uncertainties on the measurement of guH can be expressed as follows: 
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The default value for L assumed through the whole study is 1000 fh ^. This large 
value is quite essential to maintain the statistical uncertainty at the level of a few per-cent. 

\^f is the relative uncertainty on the residual background normalisation. It is mostly due 

to badly known differential cross-sections in weakly populated phase space corners. It is 
sizeable and moreover difficult to estimate. However, it pairs will be copiously produced 
at the LC, allowing a complete characterisation of them. For example, studying these 
pairs during a run at an energy where the itH process is negligible will allow to improve 
their simulation by event generators. Therefore, this uncertainty should not exceed 10%. 
We will thus repeat the analysis for two values of this uncertainty, 10% and 5%. In the 
systematic uncertainty, we just take into account the one which arises from the effective 
background normalisation since it is by far the largest one among those we can estimate 
now. 

The sensitivity factors Sstat and Ssyst in relations 3.1 and 3.2 express the dependence 
of the cross-section on the coupling squared: 
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As pointed out in section 2, the contribution from Higgs radiation off the Z to the 
signal cross-section is very small. In order to calculate the sensitivity factors, we will thus 
neglect it, allowing a very simple calculation. In this approximation, we can write: 



^tlH -gttHF'{MH,mt,s) 



And thus: 
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where s is the squared coUision energy. 
The sensitivity factors eventually read: 
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The values of Sstat and Ssyst are summarised in table 2. In this approximation, Ssyst 
is independent of the Higgs mass and of the center-of-mass energy. The results obtained 
for Mh = 120 GeV/c^ agree with those presented in [4]. 



Mh (GeV/c2) 


Sstat (fb^/^) 


^syst 


120 


3.13 


1.98 


130 


2.92 


1.98 


140 


2.72 


1.98 


150 


2.54 


1.98 


170 


2.21 


1.98 


200 


1.77 


1.98 



Table 2: Sensitivity factors for various Higgs mass values at ^/s = 800 GeV. 



4 Analysis and simulation 

4.1 Background 

The resonant background processes considered in the analysis are listed in table 3, together 
with their cross-sections (as given by the CompHEP software). 



Elementary process 


a(fb) 


e^e~ -^ qq 


1557.7 


e+e~ -^ tt 


297.3 


e+e- -^ W+W- 


4298 


e+e- -^ ZZ 


239.8 


e+e- -^ ttZ 


4.3 



Table 3: Resonant background cross-sections for ^/s = 800 GeV from CompHEP. Initial 
state radiation and beamstrahlung are taken into account. 

The processes e+e- -^ qq and e+e- -^ W+W^ exhibit a very different topology from 
the signal, however, their huge cross-section (two or three orders of magnitude larger than 
the cross-section of the signal) forbids to neglect them. 



Although its cross section is much lower, the above statements hold for the process 
e+e" -^ ZZ. 

The process e^e~ -^ iiZ and its cross section are very close to those of the signal. 

Finally, the process e^e~ -^ it has a rather large cross-section and will often mimic 
the signal. It is expected to be the main background. 

The backgrounds listed in table 3 don't account for 6 fermion production, which may 
turn out to be not negligible. The generation of these events is however a complicated and 
very time consuming task, which leads to rather large uncertainties. As we do not want 
to completely neglect these final states, we adopt the following procedure, which acts as a 
compromise. The selection is first optimised for the suppression of resonant background 
(which is by far the main background). The residual contamination from the dominant 6 
fermion processes is then estimated with the same selection criteria. The results will be 
presented with and without inclusion of the 6 fermion processes in the overall background. 

4.2 Generation of events and simulation of the detector 

The itH and iiZ partonic events were generated with CompHEP V.41.10. This program 
allows to include the initial state radiation and the beamstrahlung. These events were 
then treated with PYTHIA V.6.158 [8] for hadronization, decay and final state radiation. 

Resonant backgrounds were generated with PYTHIA V.6.158. In this case, the initial 
state radiation was considered in the structure function approach and the beamstrahlung 
was implemented with CIRCE [9]. 

The 6 fermion process partonic events were generated with WHIZARD V.1.2x^ [10]. 
As for iiH and iiZ events, they were treated with PYTHIA V.6.158 for hadronization, 
decay and final state radiation. Beamstrahlung and initial state radiation are handled by 
WHIZARD. 

Some of these 6 fermion processes receive contributions from the reactions listed in 
table 3. For these processes, the corresponding diagrams are removed from the calculation 
to avoid double counting, keeping however the interference between the resonant and non- 
resonant diagrams. This short-cut is motivated by the complexity a more correct treatment 
would introduce; it is justified by the magnitude of the interferences considered. 

The cross-sections of the 6 fermion processes are in general calculated with quite good 
accuracy. However, when coming to the event generation, things get substantially more 
complicated and time consuming. Due to CPU limitations, the number of events in the 6 
fermion samples is sometimes modest and, conversely, the differential cross sections have 
sizeable statistical uncertainties. Furthermore, some ambiguities may arise when pairing 
the particles in the final state (needed for hadronization and final state radiation) e.g. 
when there are 2 pairs of identical particles. 

The contaminations by the 6 fermion processes will thus not be estimated as precisely 
as the ones by the resonant backgrounds. The loss of resolution arising from these events 



^Versions 1.22 to 1.24 were used. 



should thus be taken as a rough estimate. This shoud not matter as the residual 6 fermion 
background is small. 

The events were further processed by SIMDET V.4 [11], the fast simulation program of 
the TESLA [12] detector, in order to take into account detector and event reconstruction 
effects. 

4.3 b-tagging 

The b-tagging is an essential tool for this analysis, especially when assuming that the 
Higgs boson decays into a pair of b-quarks. b and c-tagging algorithms developed for LEP 
and SLC experiments were adapted to TESLA and their performances studied [13] [14]. 
Recently, these tools were made available with the fast simulation of the detector [15]. 
Various algorithms are combined, including the SLD-ZVTOP vertex finder, in a neural 
network. One of the outputs is the b-probability of a jet. 

4.4 Definition of the variables used in the analysis 
4.4.1 Fox- Wolfram moments 

The Fox- Wolfram moments Hi, I = 0,1,2, .. ., are defined by 



Hi = ^^-^^Pi{cose,,), (4.1) 



e: 

where 9ij is the opening angle between hadrons i and j and -Evis the total visible energy 
of the event. The Pi{x) are the Legendre polynomials. In this paper, the moments will be 
normalized to Hq, i.e. Hiq = Hi/Hq. 

4.4.2 Heavy and light jet masses 

The particles of an event are divided into 2 classes and the invariant mass of each class 
is calculated. We obtain the heavy and light jet masses when the assignment of particles 
to the classes is such that the quadratic sum of the invariant mass of the 2 classes is 
minimised. In events with resonances (e.g. e+e^ -^ W^W~), these quantities tend to 
peak at their invariant mass. 

4.4.3 Other variables 

E-''^ . . s is the energy of the most (least) energetic iet of an event. The minimum iet 

max(mm) '^^ \ / o j j 

multiplicity is the number of objects (charged as well as neutral) in the jet which has the 
smallest multiplicity. P^'^ (i) is the probability of the i^^ jet to be a b-jet, the P^'^ (i) being 
sorted out in decreasing order, i.e. from the most b-like to the least b-like. 
The other variables are self-explanatory. 

10 



4.5 Neural network analysis 

A neural network is used to optimize the selection of events. We used the MLPfit program 
[16], a multi-layer perceptron with error backpropagation. 

5 Study of the H ^ bb decay mode 

5.1 Introduction 

When the Higgs boson decays into 66 pairs, 3 classes of final states occur. Among them, 
two can potentially allow to measure the Yukawa coupling well: the semileptonic final 
state {tt -^ W^bW^b -^ 2b2qlv with BR{tt -^ 2b2qlv) k, 43.9%) and the hadronic final 
state {it -^ W+bW-b -^ 264g with BR{tt -^ 264g) « 45.6%). These two channels are 
characterized by a large particle and jet multiplicity and an isotropic topology. The pres- 
ence of four b-jets will allow the construction of very discriminating variables. However, 
the event rate is really tiny in comparison with the background and the very crowded 
environnement will degrade clustering and b-tagging algorithms and will make invariant 
mass constraints less effective. Furthermore, hard gluon radiation combined with gluon 
splitting to 66 will often allow the background events to fake the signal. 

For both channels, a preselection sequential procedure is first applied in order to re- 
move most of the background, while keeping a high selection efficiency for the signal. A 
neural network analysis will then be performed in order to optimally use the information 
contained in the distributions of the final state characteristics. The analysis will be re- 
peated for the following Higgs boson masses: 120 GeV/c^, 130 GeV/c^, 140 GeV/c^ and 
150 GeV/c^. As the distributions hardly change in this mass window, the same procedure 
will be applied for each value of the Higgs boson mass. 

5.2 Semileptonic Channel 
5.2.1 Introduction 

The final state follows from the process: 

e+e" -^ tiH -^ W+bW~bbb -^ Ab2qlu. 

This channel is thus characterized by 4 b-jets, 2 light quark jets, one prompt lepton and 
missing 4-momentum. It has a little less statistics than the hadronic channel but the final 
state is cleaner and the presence of an isolated lepton together with missing 4-momentum 
will allow to construct powerful selection variables. Notice that, here, the hadronic channel 
{iiH -^ 464(7), the fully leptonic channel {iiH -^ Ab2l2u) and the channels where the Higgs 
boson doesn't decay into b quarks act as background processes. 
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5.2.2 Sequential analysis 

First, we require the presence in the event of at least one charged lepton (a /i^ or a e^). 
Then, we request (figure 5 and 6^): 

• 500 GeV/c2< Total visible mass < 750 GeV/c^ 

• Total multiplicity > 110 

• Number of jets (including possible isolated leptons; JADE algorithm with fjcut = 
1.10-3) > 7 

• Thrust < 0.85 

• Light jet mass > 50 GeV/c^ 

• Heavy jet mass > 150 GeV/c^ 

• Fox- Wolfram moment hlO < 0.2 

• Fox- Wolfram moment h20 < 0.6 

• Fox- Wolfram moment h30 < 0.4 

• Fox- Wolfram moment h40 < 0.5 

Now, an energetic and isolated charged lepton has to be identified. Among all the 
charged leptons (^u^ and e^) reconstructed by the detector, the one which maximises : 
El * {\ — cos9), El being the energy of the lepton and 9 the angle between its direction 
and that of the closest jet when we force the rest of particles in a 6 jet configuration (with 
JADE), is chosen. This procedure allows to choose a lepton from a leptonic decay of a W 
rather than from a B (or D) -meson decay which tends to be less energetic and isolated. 

Once the lepton is tagged, the remaining particles are forced into 6 jets with the JADE 
algorithm. We then flavour-tag the jets and finally we require (figure 6): 

• Mininum jet multiplicity > 3 

The overall preselection efficiencies and corresponding effective cross-sections^ are 
shown in table 4. Apart from hadronic tiH and ttH(H -^ bb) events, the backgrounds 
with the highest preselection efficiencies are, as expected, those of the ti and tiZ produc- 
tions. However, due to its relatively high cross-section, the main background after this 
preselection is by far the one due to the e~^e~ -^ it process. 



^The figures 4 to 22 are placed after the references. 

^The effective cross-section of a process X is defined as the remaining cross-section of this process after 

tot 
X 



a selection procedure: a'j/-' — a^f ' * ex; where ex is the selection efficiency for the process X 
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Final state 


Mh (GeV/c^) 


^preset (%) 


^effifb) 


tiH -^ Ab2qlu 


120 


61.8 


4.56-10-1 


tiH -^ Ah2qlu 


130 


62.6 


3.04-10-1 


tiH -^ Ah2qlu 


140 


64.0 


1.75-10-1 


tiH -^ Ab2qlv 


150 


65.3 


7.57-10-2 


ti 


- 


9.02 


26.8 


tiZ 


- 


23.5 


1.01 


WW 


- 


8.30 • 10-3 


3.56-10-1 


ZZ 


- 


7.27 • 10-2 


1.74 


10-1 


qq 


- 


6.25 • 10-3 


9.77 


10-2 


ttH -^ AbAq 


120 


40.3 


3.08 


10-1 


tiH -^ AbAq 


130 


40.6 


2.05 


10-1 


tiH -^ AbAq 


140 


40.2 


1.14 


10-1 


tin -^ AbAq 


150 


40.4 


4.86 


10-2 


ttH -^ 4b2l2u 


120 


7.28 


1.34 


10-2 


tiH -^ Ab2l2u 


130 


8.21 


9.54 


10-3 


tiH -^ Ab2l2u 


140 


9.44 


6.18 


10-3 


tiH -^ Ab2l2u 


150 


10.7 


2.97 


10-3 


tiH{H ^ 26) 


120 


30.3 


2.48 


10-1 


tiH{H ^ 26) 


130 


31.9 


3.36 


10-1 


tiH{H ^ 26) 


140 


33.3 


4.21 


10-1 


ttH{H ^ 26) 


150 


33.5 


4.65 


10-1 


Total background 


120 


4.52 • 10"^ 


28.9 


Total background 


130 


4.52 • 10-^ 


28.9 


Total background 


140 


4.51 • 10-1 


28.9 


Total background 


150 


4.51 • 10-1 


28.9 



Table 4: The H ^ bb semileptonic channel: preselection efficiencies (epresei) o,nd corre- 
sponding effective cross-sections (cr^ff) for various Higgs boson masses. The 4 top lines 
stand for the selected signal final state, while the next lines break down the different back- 
ground components (including those due to tiH events with another Higgs decay channel). 



5.2.3 Neural network analysis 

The distributions used for the preselection procedure still contain unexploited information. 
Some of them are then used to train a neural network together with other variables related 
to b-tagging, to the charged lepton and including the missing momentum. 
The complete list of variables is (figure 7 and 8) : 

• The total visible mass. 
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• The number of jets (including possible isolated leptons; JADE algorithm with i/cut = 
1.10-3). 

• The thrust. 

• The aplanarity. 

• The second Fox- Wolfram moment h20. 

• Etix - E^^^ ■ 

• The energy of the tagged lepton. 

• The invariant mass of the system made of the lepton and the missing momentum. 

• The cosine of the angle between the tagged lepton and the closest jet directions. 

Apart from the 5 first ones, all these variables were calculated once the event was 
forced to the lepton tagged plus 6 jet configuration. 

Once the neural network is trained, its weights are optimized for the separation of 
signal from background events. The distribution of the neural network output for these 
2 classes of events is shown in figure 9 for the case where Mh = 120 GeV/c^. From this 
figure, we observe that the separation is quite effective. However, due to the tiny cross- 
section of the signal process, the purity of the selected sample and the selection efficiency 
of the signal will be small and, eventually, the accuracy on the measurement of guH will 
be limited accordingly (see relations 3.1 and 3.2). 

5.2.4 Results 

The next step consists to apply a cut on the output of the neural network to further 
separate the signal and the background. The cut value is chosen such that it minimises the 
quadratic sum of the statistical and systematic uncertainties. The value of this optimal 
cut depends on the assumption made for the value of the uncertainty on the residual 
background normalisation and on the number of signal events, which in turn depends on 
the Higgs boson mass. 

The evolution of the selection efficiency of the signal e^gf"" , the purity of the selected 
sample /O*"™^'^, their products, the statistical, systematic and total uncertainties as a 
function of the cut on the neural network output are shown in figure 3 for the case where 

Mh = 120 GeV/c and — ^ff- = 5%. As we increase the cut, the purity gets higher and its 



'BG 



derivative is increasing. The systematic uncertainty behaves the same way as 
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Figure 3: The H ^ bb semileptonic channel: (left) selection efficiency, purity and their 
product, (right) statistical, systematic and total uncertainties on the measurement of quh 

as functions of the value of the cut on the neural network output for \^f = 5% and 

Mh = 120 GeV/c^. 

It is a decreasing function of p, which varies fast for values of p close to and less as p 
increases. The evolution of the selection efficiency is opposite to the one of the purity. 
Their product has thus a small variation. It however increases until a particular cut value 
(0.92 for this case). The statistical uncertainty will thus slowly decrease until this cut 
value. A higher cut provokes a sharp drop of the efficiency and a degradation of the 
statistical uncertainty. Eventually, the optimal cut is about 0.95 for the case under study. 
The selection efficiencies and the corresponding effective cross-sections are shown in 

table 5 for Mh = 120 GeV/c^ and — jf^ = 5%. In this case, the number of selected events 

amounts to 107 for the signal and to 192 for the background. The main background after 
selection is due to top-pair production. 

The uncertainties on guH are shown in table 6 and on the figure 4. They range from 

9.1% (11.7%) to 48.7% (65.9%) for ^^ = 5% (10%). The resolution degrades with 

'^BG 

increasing mass due to the reduction of the statistics of the signal. 
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Final state 


esel (%) 


CTe// {fh) 


tiH -^ Ab2qlu 


14.6 


1.07-10-1 


tt 


4.26 • 10-^ 


1.27-10-1 


tiZ 


9.36-10-1 


4.01-10-2 


WW 


< 3.7 -10-5 


< 1.6-10-3 


zz 


5.77-10-4 


1.38-10-3 


QQ 


< 10-4 


< 1.6-10-3 


ttH -^ AbAq 


1.82 


1.40-10-2 


tiH -^ Ab2l2u 


1.32 


2.33-10-3 


ttH{H ^ 26) 


9.6-10-1 


7.86-10-3 


Total background 


3.01 - 10-^ 


1.92-10-1 



Table 5: The H ^ bb semileptonic channel: selection efficiencies (egei) and corresponding 



effective cross-sections (cfeff) for 



Act 



e// 



77r 



5% and M// = 120 GeV/c^ . 



Mh (GeV/c2) 


^^'^i 
-P^^' 


signal 
^sel 


sample 
Psel 


( ^9ttH \ 

\ gttH ^Stat 


^ gttH 'syst 


AgttH 
gttH 


120 


5% 
10% 


14.6% 
10.7% 


36.0% 

39.7% 


8.0% 
8.8% 


4.5% 
7.7% 


9.1% 
11.7% 


130 


5% 
10% 


16.3% 

12.9% 


26.0% 

28.8% 


10.9% 
11.6% 


7.2% 
12.5% 


13.0% 
17.1% 


140 


5% 
10% 


17.4% 
11.1% 


17.1% 
20.1% 


17.5% 
20.2% 


12.2% 
20.1% 


21.3% 
28.5% 


150 


5% 
10% 


13.8% 
13.8% 


9.0% 
9.0% 


41.4% 
41.4% 


25.6% 
51.2% 


48.7% 
65.9% 



Table 6: The H ^ bb semileptonic channel: expected uncertainty on the measurement of 
gttH for various Higgs boson masses. Selection efficiency of the signal (e*g^"" ) and purity 
of the selected sample (p^e™*' ^ ) are also shown. 

5.2.5 Precision loss due to 6 fermion processes 

In the case of the final state under study, the most contaminating 6 fermion processes are 
the ones made of at least 4 b-quarks. They were studied for Mh = 120 GeV/c^. Their 
cross section and selection efficiency are shown in table 7 and the loss of resolution on the 
measurement of guH is shown in table 8. One observes that it can be neglected. 
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Final state 


<m 


^sel 


bbbbbb 


6.4 • 10^3 


5.0% 


bbbbqq {q = u,d) 


1.2-10-1 


1.0% 


bbbbss 


6.3-10-2 


1.3% 


bbbbcc 


5.1-10-2 


2.0% 


bbbbtt 


6.9 - 10-3 


14.2% 


bbqqtt (q = u, d, s) 


8.9 - 10-3 


5.2% 


bbcctt 


3.5 - 10-3 


6.6% 


bbtbdu * 


~ 1 - 10-2 


5.4% 


bbtb^u * 


~ 3 - 10-3 


14.2% 



Table 7: The H ^ bb semileptonic channel: cross section (a) and selection efficiencies 



(^sel) of the 6 fermion processes for 



Act' 



e// 



BG 

TT 

BG 



5% and Mh = 120 GeV/c^. For events 



marked with a *, some diagrams were removed from the calculation to avoid double count- 
ing. 



^RG 


^^ (without 6f) 


^(with 6f) 


5% 

10% 


9.1% 
11.7% 


9.3% 
12.1% 



Table 8: The H ^ bb semileptonic channel: expected uncertainty on the measurement of 
9ttH without and with inclusion of the 6 fermion processes in the analysis for Mh = 120 
GeV/c^. 

5.3 Hadronic channel 
5.3.1 Introduction 

The final state follows from the process: 



e"'"e 



ttH -^ W+bW-bbb -^ 4b4q. 



This channel is thus characterized by 4 b-jets and 4 light quark jets. It benefits from 
the highest branching ratio. However the environnement is still more crowded than in 
the semileptonic channel and the event reconstruction will be quite affected. Notice that, 
here, the semileptonic channel {tiH — >■ 4:b2qliy), the fully leptonic channel {iiH -^ 4:b2l2v) 
and the channels where the Higgs boson doesn't decay into b quarks act as background 
processes. 
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5.3.2 Sequential analysis 

We first require: 

• Total visible mass > 560 GeV/c^. 
Total multiplicity > 120. 



• 



• Number of jets (including possible isolated leptons; JADE algorithm with ycut = 
1.10-3) > 7. 

• Thrust < 0.85. 

The events satisfying these requirements are forced in a 8 jet configuration (using 
JADE) and a second set of criteria is then applied: 

• Minimum jet multiplicity > 1. 

• Minimum dijet invariant mass > 15 GeV/c^. 

• E^""^ > 20 GeV. 



ram 



All these distributions are shown in figure 10. The overall preselection efficiencies and 
the corresponding effective cross-sections are shown in table 9. The remarks made on the 
residual background in the semileptonic channel (see section 5.2.2) apply here as well. 

5.3.3 Neural network analysis 

We then train a neural network to use the information contained in the distributions in a 
more powerful way. The variables used are: 

• The total visible mass. 

• The number of jets (including possible isolated leptons; JADE algorithm with ycut = 
1.10-3). 

• The thrust. 

• The aplanarity. 

• The light jet mass. 

• The heavy jet mass. 

• The second Fox- Wolfram moment h20. 
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Final state 


Mh GeV/c^ 


^presel \/o) 


^effifb) 


tiH -^ AbAq 


120 


63.2 


4.85-10-1 


tiH -^ AhAq 


130 


64.0 


3.26-10-1 


tiH -^ AhAq 


140 


64.3 


1.81-10-1 


tiH -^ AbAq 


150 


64.6 


7.84-10-2 


ti 


- 


7.37 


21.9 


tiZ 


- 


22.8 


9.74-10-1 


WW 


- 


1.10-10-1 


4.74 


zz 


- 


2.56-10-1 


6.14-10-1 


qq 


- 


6.93 - 10-2 


1.08 


ttH -^ Ab2qlu 


120 


14.1 


1.04-10-1 


tiH -^ Ah2qlu 


130 


14.7 


7.20-10-2 


tin -^ Ab2qlu 


140 


16.0 


4.32-10-2 


tiH -^ Ah2qlu 


150 


16.4 


1.91-10-2 


tin -^ Ah2l2u 


120 


3.60-10-1 


6.66-10-^ 


tiH -^ Ab2l2u 


130 


3.60-10-1 


4.42-10-4 


tiH -^ Ab2l2u 


140 


5.20-10-1 


3.52-10-4 


tiH -^ Ab2l2u 


150 


6.85-10-1 


2.00-10-4 


tiH{H -^ 2b) 


120 


26.4 


2.16-10-1 


tiH{H -^ 2b) 


130 


29.5 


3.14-10-1 


tiH{H -^ 2b) 


140 


30.6 


3.85-10-1 


ttH{H -^ 2b) 


150 


31.9 


4.39-10-1 


Total background 


120 


4.46 - 10-1 


28.6 


Total background 


130 


4.47 - 10-1 


28.6 


Total background 


140 


4.48 - 10-1 


28.7 


Total background 


150 


4.49 - 10-1 


28.7 



Table 9: The H ^ bb hadronic channel: preselection efficiencies (epresei) CL^d correspond- 
ing effective cross- sections (cyeff) for various Higgs boson masses. The 4 top lines stand 
for the selected signal final state, while the next lines break down the different background 
components (including those due to tiH events with another Higgs decay channel) 



JZ/rt 



et 
max 



Ei 



jet 
min' 



E^ 



et 



max{min) 



and the Pi 



.jet I 



are calculated once the event is forced into 8 jets. The 



distributions of these variables (after the preselection) are shown in figures 11 and 12. 



5.3.4 Results 

The result of the training is shown in figure 13. We apply the same procedure as for 
the semileptonic channel (see section 5.2.4). An example of the values of the selection 
efficiencies and of the corresponding effective cross-sections is shown in table 10. As 
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for the semileptonic channel, the main background is due to top pair production. The 
resolutions on gun are shown in table 11 and on figure 4. They range from 8.3% (10.1%) 

to 42.7% (51.2%) for _^ff = 5% (10%). As in the case of the semileptonic channel, 



eff 



the resolution degrades with increasing mass due to the reduction of the statistics of the 
signal. 



Final state 


esel (%) 


^effifb) 


tiH -^ AbAq 


14.8 


1.14-10-1 


tt 


3.10-10-^ 


9.29-10-2 


tiZ 


1.11 


4.74-10-2 


WW 


< 3.7 -10-5 


< 1.6-10-3 


zz 


8.71 • 10-4 


2.08-10-3 


QQ 


<io-4 


< 1.6-10-3 


ttH -^ Ab2qlu 


2.44 


1.80-10-2 


tiH -^ 4b2l2iy 


3.14-10-2 


5.55-10-5 


ttH{H ^ 26) 


7.9-10-1 


6.47-10-3 


Total background 


2.61 - 10-'^ 


1.67-10-1 



Table 10: The H ^ bb hadronic channel: selection efficiencies (egei) o-nd corresponding 



effective cross- sections (cfeff) for 



Act 



BG 



5% andMH = 120 GeV/c^ 



Mh (GeV/c2) 


^ 

'^p" 


signal 
^sel 


sample 
Psel 


( ^9ttH \ 

\ gttH 'Stat 


^ gttH 'syst 


^gttH 
gttH 


120 


5% 
10% 


14.8% 
12.0% 


40.5% 

44.7% 


7.4% 

7.9% 


3.7% 
6.3% 


8.3% 
10.1% 


130 


5% 
10% 


15.6% 

9.5% 


27.5% 
35.5% 


10.8% 

12.2% 


6.7% 
9.2% 


12.7% 
15.3% 


140 


5% 
10% 


13.7% 
8.9% 


19.9% 

23.7% 


18.3% 
20.7% 


10.2% 
16.3% 


20.9% 
26.4% 


150 


5% 
10% 


12.4% 
8.8% 


11.7% 
14.2% 


38.2% 
41.1% 


19.1% 
30.5% 


42.7% 
51.2% 



Table 11: The H ^ bb hadronic channel: expected uncertainty on the measurement of 
9ttH for various Higgs boson masses. Selection efficiency of the signal (e*g^"" ) and purity 
of the selected sample (p^g["^ ^J are also shown. 
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5.3.5 Precision loss due to 6 fermion processes 

The 6 fermion processes which may degrade the resolution on guH in this channel are 
the same as in the H ^ bb semileptonic channel and they were also studied for Mh = 
120 GeV/c^. Their cross sections and selection efficiencies are shown in table 12. The 
precisions obtained on gun when the 6 fermion background is included are shown in 
table 13. As in the semileptonic channel, the degradation of the measurement due to the 
6 fermion background is negligible. 



Final state 


<m 


fse/ 


bbbbbb 


6.4 • 10-3 


5.8% 


bbbbqq (q = u,d) 


1.2-10-1 


1.6% 


bbbbss 


6.3-10-2 


1.8% 


bbbbcc 


5.1-10-2 


2.0% 


bbbbtt 


6.9 - 10-3 


31.8% 


bbqqti (q = u, d, s) 


8.9 - 10-3 


12.5% 


bbcctt 


3.5 - 10-3 


17.2% 


bbtbdu * 


~ 1 - 10-2 


14.0% 


bbtbfiV * 


~ 3 - 10-3 


1.2% 



Table 12: The H ^ bb hadronic channel: cross section (a) and selection efficiencies (csei) 



of the 6 fermion processes for 









5% and Mjj = 120 GeV/(P. For events marked 



with a *, some diagrams were removed from the calculation to avoid double counting. 





^^ (without 6f) 


^^"^ (with 6f ) 


5% 
10% 


8.3% 
10.1% 


8.5% 
10.5% 



Table 13: The H ^ bb hadronic channel: expected uncertainty on the measurement of 
guH without and with inclusion of the 6 fermion processes in the analysis for Mh = 120 
GeV/c^. 



6 Study of the H -^ W^W decay mode 

6.1 Introduction 

As we can see from the previous section, the measurement of gun in the H ^ bb decay 
mode degrades quite rapidly as the Higgs boson mass increases. The reasons are the drop 
of the cross-section and even more, the decrease of the branching ratio of the Higgs boson 
into bb for the benefit of the H -^ W^W~ decay. We will thus try to exploit this latter 
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mode when the Higgs boson is heavier than 140 GeV/c^, and see whether it also improves 
the accuracy on quh for hghter Higgs boson masses. 

In this decay mode, four intermediate W bosons are present, leading to several classes 
of final states. Unlike the final states where the Higgs boson decays into a pair of b-quarks, 
there are only two b-jets in the event, thus the b-tagging is no longer an essential point of 
the analysis. Therefore, particular final states have to be identified which can allow good 
signal and background separation and which have enough statistics. Two such channels 
were found: the "2 like sign lepton plus 6 jet channel", when two W bosons of the same 
sign decay leptonically whereas the two remaining ones decay hadronically, and the "single 
lepton plus 8 jet channel" when only one of the W's decays leptonically. As for channels 
where the Higgs boson decays into pairs of b-quarks, these final states have large particle 
and jet multiplicities, an isotropic topology and a tiny event rate. For each channel, the 
same analysis will be repeated for 4 values of the Higgs boson mass within the range 130 
GeV/c2 - 200 GeV/c2. 

6.2 The 2 like sign lepton plus 6 jet channel 

6.2.1 Introduction 

The final state follows from the process: 

e+e" -^ tiH -^ W+bW-bW+W -^ 2l^2v2bAq 

Its branching ratio is: 

BR{ttH ^2l^2v2bAq) = BR{H ^W+W~) 

*BR{W^W^ -^ Aq) * BR{W^W^ -^ 2l2v) * 2 
^ 9.Q%*BR{H ^W^W) 

This channel is thus characterized by a missing 4-momentum, two energetic and iso- 
lated charged leptons of the same sign, 4 light quark jets and 2 b-jets. In comparison 
with the case where no restriction on the charged lepton signs is made, requiring the two 
leptons to have the same sign divides the branching ratio of the signal by a factor 3 but 
the background can be very effectively suppressed. 

Here, ttH{H -^ WW) -^ 2l^2v2bAq and ttH{H -p> WW) act as background processes. 
The analysis will be purely sequential, no neural network being applied. 

6.2.2 Analysis 

We first apply a set of criteria related to topological variables (figure 14): 

• 400 GeV/c^ < Total visible mass < 700 GeV/c^. 

• 85 < Total multiphcity < 160. 
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• Number of jets (including possible isolated leptons; JADE algorithm with ycut = 
1.10-3) > 6. 

• Light jet mass > 100 GeV/c^. 

• Heavy jet mass > 150 GeV/c^. 

• Fox- Wolfram moment hlO < 0.2. 

• Fox- Wolfram moment h20 < 0.45. 

• Fox- Wolfram moment h30 < 0.3. 

• Fox- Wolfram moment h40 < 0.3. 

The particles (discarding the charged leptons) of the surviving events are forced to a 
6 jet configuration (using JADE) and the jets are flavour-tagged. The events are then 
required to fulfill the following criteria (figure 14): 

• Minimum dijet invariant mass > 15 GeV/c^. 

• Pt\l) > 0.2. 

The request of only one b-jet allows to preserve a high selection efficiency. Now, we 
make use of the lepton content of the signal events to further eliminate the background. 
For each charged lepton reconstructed in the event, we calculate the transverse momentum 
with respect to the closest jet. Next, we classify the leptons from the most isolated to the 
least isolated according to this quantity. In the event, there should be 2 and only 2 isolated 
leptons. Moreover, they should have the same sign. As they come from the primary vertex, 
the significance of their impact parameter^ should be small in comparison with the one of 
a lepton coming from a b- or c-meson decay. Therefore, we require (figure 15): 

• The 2 most isolated leptons have the same sign. 

• A pt between the most isolated lepton and the other jets > 5 GeV/c. 

• A pt between the second most isolated lepton and the other jets > 5 GeV/c. 

• A Pt between the third most isolated lepton and the other jets < 5 GeV/c. 

• -0.001 < Significance of most isolated lepton impact parameter < 0.001. 

• -0.001 < Significance of second most isolated lepton impact parameter < 0.001. 



*For a charged track, the impact parameter is the distance of closest approach between this track and 
the primary vertex and the significance of the impact parameter is the ratio between the impact parameter 
and its estimated uncertainty. 
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The overall selection efficiencies and the corresponding effective cross-sections are 
shown in table 14. The main background is again due to top pair production. As the 
number of signal events is small in this channel but the purity high, the total uncertainty 
on guH is dominated by the statistical uncertainty. Therefore, the relative uncertainty 
on the residual background normalisation has only a modest influence on the total uncer- 
tainty. As a matter of fact, the procedure is not reoptimised for each value of the relative 
uncertainty on the residual background normalisation. 



Final state 


esel (%) 


'^effifb) 


tiH -^ 2l^2v2h4.q; Mr = 130 GeV/c^ 


20.2 


1.26 -10-2 


tiH -^ 2l^2u2hAq; Mh = 150 GeV/c^ 


24.6 


2.72 -10-2 


iiH -^ 2l^2u2bAq; Mh = 170 GeV/c^ 


25.5 


2.96 -10-2 


tiH -^ 2l^2u2b4q; Mh = 200 GeV/c^ 


26.3 


1.47 -10-2 


tt 


4.40 • 10"^ 


1.31-10-2 


tiZ 


7.33-10-2 


3.14-10-3 


WW 


7.41 • 10-^ 


3.18-10-3 


ZZ 


1.02 • lO-''^ 


2.44-10-3 


qq 


< 8.33 • 10-5 


<1.30-10-3 


ttH{H -^ WW) ^ 2l^2v2bAq; Mh = 150 GeV/c^ 


2.65 • 10-1 


2.77-10-3 


ttH{H -^ WW); Mh = 150 GeV/c^ 


8.0-10-2 


3.94-10-4 


Total background ; Mh = 150 GeV/c^ 


3.84 • 10-4 


2.46-10-2 



Table 14: The two like sign lepton plus 6 jet channel: selection efficiencies (egei) and 
corresponding effective cross-sections ((Jeff) for various Higgs boson masses. The 4 top 
lines stand for the selected signal final state, while the next lines break down the differ- 
ent background components (including those due to tiH events with another Higgs decay 
channel). 



6.2.3 Results 

The resolutions on the measurement of gun are shown in table 15 and on figure 4. They 



range from 12.8% (13.4%) to 24.9% (26.3%) for 



e// 



eff 



5% (10%). For the lowest values 



of the Higgs mass range, the resolution is not good as the branching ratio of the Higgs into 
pairs of W's is very small. For Mh ~ 170 GeV/c2, the resolution curve has a minimum 
which corresponds to the maximum of the branching ratio. Then, for higher mass values, 
the resolution degrades due to the decrease of the cross section and of the branching ratio. 



6.2.4 Precision loss due to 6 fermion processes 

The 6 fermion processes which may degrade the resolution on gun in this channel are 
listed in table 16 with their cross section and selection efficiency for Mh = 150 GeV/c2. 
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Mh (GeV/c2) 




signal 
^sel 


sample 
Psel 


( AgttH \ 


'^ gttH '^yst 


AgttH 
gttH 


130 


5% 
10% 


20.2% 
20.2% 


33.9% 
33.9% 


24.4% 
24.4% 


4.9% 
9.9% 


24.9% 
26.3% 


150 


5% 
10% 


24.6% 
24.6% 


52.1% 
52.1% 


13.4% 
13.4% 


2.3% 
4.6% 


13.6% 
14.2% 


170 


5% 
10% 


25.5% 
25.5% 


54.0% 
54.0% 


12.7% 
12.7% 


2.2% 
4.3% 


12.8% 
13.4% 


200 


5% 
10% 


26.3% 
26.3% 


37.3% 
37.3% 


21.6% 
21.6% 


4.3% 
8.5% 


22.0% 

23.2% 



Table 15: The two like sign lepton plus 6 jet channel: expected uncertainty on the mea- 
surement of gttH for various Higgs boson masses. Selection efficiency of the signal (e^g^"" ) 
and purity of the selected sample (Pgg™^ '^ ) are also shown. 

The loss of resolution on gttH when the 6 fermion background is included in the analysis 
is shown in table 17, showing that it may be neglected. 



Final state 


<fb) 


fse/ 


bbude^Vf, * 


~ 1.5 


8.4 • 10-3% 


bbtbdu * 
bbtbfiV * 


~ 6 • lO"-'^ 
~ 1.5 -10-3 


< 0.2% 
0.4% 


tiqqq'q' {q, q' = u, d, s) 

qqtblVi {q = u, d, s,c; I = e~, /i^) * 

qqtbdu {q = u,d,s,c) * 


1.2-10-^ 
~ 7.6 -10-2 

~ 1.3 -10-1 


0.2% 
0.2% 
0.05% 



Table 16: The two like sign lepton plus 6 jet channel: cross section (a) and selection 
efficiencies (egei) of the 6 fermion processes for — j^ = 5% and Mh = 150 GeV/c . For 






events marked with a *, some diagrams were removed from the calculation to avoid double 
counting. 





^^"« (without 6f ) 
gttH ^ ' 


^^"^ (with 6f ) 
gttH ^ ' 


5% 
10% 


13.6% 

14.2% 


13.9% 
14.5% 



Table 17: The two like sign lepton plus 6 jet channel: expected uncertainty on the mea- 
surement of gttH without and with inclusion of the 6 fermion processes in the analysis for 
Mi/ = 150 GeV/c^. 
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6.3 The single lepton plus 8 jet channel 

6.3.1 Introduction 

The final state follows from the process: 

e+e- -^ tiH -^ W+bW-bW+W -^ lu2b&q 
Its branching ratio is: 

BR{ttH ^ lv2bQq) = 4.* BR{H ^W^W) * BR{W ^ Iv) * {BR{W ^ 2q)f 
« AQ7o*BR{H^W+W-) 

This channel is thus characterized by a missing 4-momentum, one prompt charged 
lepton, 6 light quark jets and 2 b-jets. This signature is less singular than the one of the 
previous channel but the branching ratio is about 4 times larger. This final state is close to 
the one of the H ^ bb semileptonic channel, the analysis will thus be very similar. Notice 
that, here, ttH{H -^ WW) ^^ lv2bQq and ttH{H ^^ WW) act as background processes. 

6.3.2 Sequential analysis 

We first request the presence in the event of at least one charged lepton (a /i^ or a e^). 
Then, we require (figure 16 and 17): 

• 500 GeV/c2< Total visible mass < 750 GeV/c^. 

• Total multiphcity > 110. 

• Number of jets (including possible isolated leptons; JADE algorithm with ycut = 
1.10-3) > 8. 

• Thrust < 0.8 

• Fox- Wolfram moment hlO < 0.1. 

• Fox- Wolfram moment h20 < 0.5. 

• Fox- Wolfram moment h30 < 0.3. 

• Fox- Wolfram moment h40 < 0.3. 

• Light jet mass > 100 GeV/c^. 

• Heavy jet mass > 150 GeV/c^. 
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To identify an energetic and isolated charged lepton, we proceed as in the H ^ bb 
semileptonic channel. The only difference is that, here, to tag the lepton as well as when 
it has been tagged, the remaining particles are forced to a 8 jet configuration (with the 
JADE algorithm). Then, the jets are flavour-tagged and finally, we require (figure 17): 

• Mininum jet multiplicity > 3. 

The overall preselection efficiencies and the corresponding effective cross-sections are 
shown in table 18. The main background after this preselection comes from the ti pro- 
duction. 



Final state 


Cprese/ (%) 


^e//(/fc) 


tiH -^ Iiy2b6q; Mh = 130 GeV/c^ 


50.3 


1.31-10-^ 


tlH -^ Iiy2b6q; Mh = 150 GeV/c^ 


55.6 


2.57-10-1 


tiH -^ Iu2b6q; Mh = 170 GeV/c^ 


58.3 


2.83-10-1 


tlH -^ lu2bQq- Mh = 200 GeV/c^ 


60.8 


1.42-10-1 


tt 


2.76 


8.22 


ttz 


13.8 


5.91-10-1 


WW 


1.73-10^2 


7.45-10-1 


zz 


6.55-10^2 


1.60-10-1 


qq 


2.33 • 10-3 


3.63-10-2 


ttH{H -^ WW) -^ lv2bQq; Mh = 150 GeV/c^ 


18.9 


1.31-10-1 


ttH{H ^ WW); Mh = 150 GeV/c^ 


26.5 


1.30-10-1 


Total background; Mh = 150 GeV/c^ 


1.54 -10^^ 


10.0 



Table 18: The single lepton plus 8 jet channel: preselection efficiencies (epresel) '^''^d cor- 
responding effective cross- sections (c^eff) for various Higgs boson masses. The 4 top lines 
stand for the selected signal final state, while the next lines break down the different back- 
ground components (including those due to tiH events with another Higgs decay channel) 



6.3.3 Neural network analysis 

The distributions of some of the preselection variables still contain unused information. 
They are recycled to train a neural network together with other variables related to b- 
tagging, to the prompt lepton and to the missing momentum. 

The complete list of variables used in the network is (figure 18 and 19): 

• The total visible mass. 

• The number of jets (including possible isolated leptons; JADE algorithm with i/cut = 
1.10-3). 
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The thrust. 

The aplanarity. 

The energy of the tagged lepton. 

The invariant mass of the system made of the lepton and the missing momentum. 

The cosine of the angle between the tagged lepton and the closest jet directions. 



E; 



,jet 



E: 



,jet 



max min ' 

• Maximum dijet invariant mass— minimum dijet invariant mass. 

• The second Fox- Wolfram moment h20. 

• The light jet mass. 

• The heavy jet mass. 

The lepton related variables are calculated once the event has been forced to the 
1 lepton plus 8 jet configuration, as well as the variables E^ax, -^miw maximum dijet 
invariant mass and minimum dijet invariant mass. 

6.3.4 Results 



Final state 


esel (%) 


^effifh) 


ttH -^ Iu2b6q 


16.7 


7.71-10-^ 


it 


5.65-10^2 


1.68-10-1 


ttz 


4.77-10-1 


2.04-10-2 


WW 


1.48 • 10"'' 


6.37-10-3 


zz 


3.02 • 10"^ 


7.25-10-3 


qq 


< 10-^ 


< 1.6 - 10-3 


ttH{H -^ WW) -^ Iu2b6q 


1.24 


8.59-10-3 


ttH{H -^ WW) 


2.35 


1.16-10-2 


Total background 


3.30 • 10-^ 


2.22-10-1 



Table 19: The single lepton plus 8 jet channel: selection efficiencies (tsei) cmd correspond- 



ing effective cross-sections (cTeff) for 



eff 






5% and Mh = 150 GeV/(?. 



The result of the training session is shown in figure 20. We apply the same procedure as 
for the H ^ bb semileptonic channel (see section 5.2.4) and an example of the values of the 
selection efficiencies and the corresponding effective cross-sections are shown in table 19. 
The main background after selection is again due to top-pair production. The precisions 
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we get for guH are shown in table 20 and on figure 4. They range from 10.2% (12.1%) to 



35.2% (47.2%) for 



Act' 



e// 



eff 



5% (10%). The remarks made for the previous channel apply 



here as well. For the lowest values of the Higgs mass range, the resolution is not good as 
the branching ratio of the Higgs into pairs of W's is very small. For Mh ~ 170 GeV/c^, 
the resolution curve has a minimum which corresponds to the maximum of the branching 
ratio. Then, for higher mass values, the resolution degrades due to the decrease of the 
cross section and of the branching ratio. 



Mh (GeV/c2) 




signal 
^sel 


sample 
Psel 


( ^9ttH \ 

\ gttH 'stat 


( ^9ttH \ 

^ 9ttH 'syst 


AgttH 

9ttH 


130 


5% 
10% 


8.9% 

8.9% 


12.2% 
12.2% 


30.1% 

30.1% 


18.2% 
36.4% 


35.2% 

47.2% 


150 


5% 
10% 


16.7% 
9.3% 


25.8% 
33.4% 


11.3% 
13.4% 


7.3% 
10.1% 


13.5% 
16.8% 


170 


5% 
10% 


15.9% 
12.5% 


38.0% 

42.5% 


9.3% 
10.0% 


4.1% 
6.8% 


10.2% 
12.1% 


200 


5% 
10% 


16.0% 
12.5% 


27.5% 
32.2% 


15.7% 
16.5% 


6.6% 
10.6% 


17.1% 
19.6% 



Table 20: The single lepton plus 8 jet channel: expected uncertainty on the measurement of 
9ttH for various Higgs boson masses. Selection efficiency of the signal (e^g^"" ) and purity 
of the selected sample (p^g["^ '^) are also shown. 



6.3.5 Precision loss due to 6 fermion processes 

The 6 fermion processes which may degrade the resolution on guH in this channel are the 
same than in the previous channel. They are listed in table 21 with their cross section 



Final state 


<^{fb) 


f^sel 


bbude^Vf, * 


~ 1.5 


0.04% 


bbtbdE * 
bbtbfiV * 


~ 6 • 10-3 
~ 1.5 -10-3 


2.4% 
3.4% 


tiqqq'q' {q, q' = u, d, s) 

qqtblV [q = u,d,s,c; I = e~,fi~) * 

qqtbdu {q = u, d, s, c) * 


1.2-10-^ 
~ 7.6 • 10-2 
~ 1.3 -10-1 


7.6% 
1.7% 
1.2% 



Table 21: The single lepton plus 8 jet channel: cross section (a) and selection efficiencies 



(esel) of the 6 fermion processes for 



Act' 



e// 



' BG 



5% and Mh = 150 GeV/c^. For events 



marked with a *, some diagrams were removed from the calculation to avoid double count- 
ing. 
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and selection efficiency for Mh = 150 GeV/c^. The loss of resolution on quh when the 
6 fermion background is included in the analysis is shown in table 22. This channel is 
slightly more affected by the 6 fermion background than the previous ones because its 
signature is less distinctive. The loss of resolution is however still very small. 



A^^.^! 

<^'^^. 


^^ (without 6f) 


^7^ (with 6f) 


5% 
10% 


13.5% 
16.8% 


~14% 
~ 17.5% 



Table 22: The single lepton plus 8 jet channel: expected uncertainty on the measurement 
of QttH without and with inclusion of the 6 fermion processes in the analysis for Mh = 150 
GeV/c^. 



7 Influence of the limited knowledge on the SM input pa- 
rameters 

The results presented in this note were obtained for fixed values of the SM input param- 
eters. The question was addressed whether the present limited knowledge on these values 
has any significant influence on the expected precision on guH- 

The largest effect, if any, is expected from the present experimental uncertainty on 
the top quark mass (which is of the order of 5 GeV/c^). A change in nit modifies sev- 
eral ingredients of the analysis: the selection efficiencies for signal and background, the 
top-pair production cross-section (cj^j), the signal cross-section (cr^j//) and the function 
F{MH,rnt,s) entering equation 3.5 and 3.6, a.s.o. 

The study presented here was performed by moving mt from its central value of 175 
GeV/c^ to 170 and 180 GeV/c^. The change in a^jj^ was observed to be largest for low 
values of Mh (i.e. 120 GeV/c^), where it amounted to it 5 %, while the main background 
cross-section (a^j) varied by it 1 % only. Signal events were generated with mt = 170 
and 180 GeV/c^, Mh being fixed to 120 GeV/c^, and the analysis of the Higgs decay into 
b-quark pairs was repeated for the so-called semileptonic channel (see section 5.2). The 
distributions of several variables entering the event selection are displayed on figure 21 and 
22 for each value of mt- The changes observed are tiny, illustrating that the preselection 
and selection efficiencies remain essentially unchanged. Overall, the study demonstrated 
that the precision on gttn varies by a negligible amount when moving mt from 170 to 180 
GeV/c^ . The conclusion is therefore that the present accuracy on the SM input parameters 
has no significant influence on the predicted accuracy on gttH- 
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8 Conclusion and outlook 

The process e"^e~ -^ tiH allows in principle a direct measurement of the top-Higgs Yukawa 
coupling. We presented a realistic feasability study in the context of a future e~^e~ Linear 
Collider such as TESLA, that takes into account all dominant physical backgrounds, the 
main radiative effects (initial and final state radiation and beamstrahlung) as well as 
detector and event reconstruction effects. A rough estimate of the loss of precision due to 
6 fermion background processes was presented. The effect of the uncertainty on the top 
quark mass was also addressed. The e~^e~ -^ e~^e~jj -^ e"''e~hadrons events, which may 
superimpose on physics events, were neglected. However, the impact of this background 
was studied for the e"'"e~ -^ W~^W~vu -^ Hi/v process and found negligible [18]. Four 
channels were studied and the analysis repeated for several values of the Higgs boson mass 
ranging from 120 GeV/c^ to 200 GeV/c^. A collision energy of 800 GeV and an integrated 
luminosity of 1000 fb^^ were assumed. The accuracies which are obtained account for the 
statistical uncertainty and for the systematic uncertainty arising from a limited knowledge 
of the background normalisation. 

For Higgs boson masses under ~ 135 GeV/c^, the main decay mode is H ^ bb. Two 
channels were analysed - the hadronic and the semileptonic channels - and we presented the 
results for the Higgs mass range: 120 GeV/c^ - 150 GeV/c^. The measurement precision 
found for a mass of 120 GeV/c^ in the H ^ bb mode is slightly worse than in [4], due to 
refinements of the present analysis which make it more realistic. The resolution degrades 
with increasing mass, due to the reduced event rate. In this decay mode, the ability to 
identify b-jets is of major importance. 

For Higgs boson masses above ~ 140 GeV/c^, the H -^ W^W^ decay mode yields 
higher precision than the H ^ bb one as the number of signal events in the latter mode 
becomes too tiny. Two channels were analysed - the 2 like sign lepton plus 6 jet channel 
and the single lepton plus 8 jet channel - and we presented the results for the Higgs mass 
range: 130 GeV/c^ - 200 GeV/c^. 

We showed that the 6 fermion background has a very modest influence on the mea- 
surement of the coupling, and that the present limited knowledge on the Standard Model 
input parameters has no significant influence on the predicted accuracy on guH- 

As a final result, the four channels studied in this paper are combined^ to get the 
global precision and these results are shown in table 23 and on figure 4. 

The expected accuracy on guH is better than « 10% over most of the mass range (up 
to Mff ~ 180 GeV/c^), even if the knowledge of the background normalisation is only 

at the 10% level. In the most favourable case {Mh= 120 GeV/c and —^f- = 5%), the 

accuracy is about 6%. It is also good (8 to 9%) for Mh « 170 GeV/c^, corresponding to 
the maximum of the branching ratio of the Higgs boson into pairs of W's. For the less 
favourable case {Mh= 200 GeV/c^), the accuracy is however still better than 15% even if 

^We combine the uncertainties found without inclusion of 6 fermion processes as the degradation due 
to them has only been roughly estimated and is anyway small. 
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'^BG 



10%. 



We observe that the best resolutions are obtained for the lowest Higgs mass values, 
which are the most likely as predicted by precision measurements. 

At the LHC, under some conditions which make the analysis more model-dependent, 
the expected accuracy on quh with 300 fb^^ lies in the range 10% — 20% for Mh G 
[100 GeV/c^-200 GeV/c^] [19]. The precise measurements of the Higgs branching ratios 
available already at a low energy run of the Linear Collider will allow to remove the model- 
dependence of the LHC results [20]. The achievable accuracies are however not as good 
as those expected at a high energy run of the Linear Collider. 

Finally, it should be stressed that there is certainly room for substantial improvement of 
the study exposed in this note. For instance, the b-tagging does not include all observables 
(e.g. vertex charge), more efficient jet reconstruction and particle flow algorithms could 
improve the reconstruction, other final states could be included, the reconstruction of 
r's could enhance the number of signal events in leptonic channels and the analysis itself 
(selection criteria, neural network inputs, training and architecture) can also be optimized. 
Moreover, the analysis may be extended above M/^ =200 GeV/c^ in order to cover the full 
Higgs mass range allowed by precision measurements. 



Mh (GeV/c2) 


'^BG 


^guH 
guH 


120 


5% 
10% 


6.1% 
7.6% 


130 


5% 
10% 


8.3% 
10.2% 


150 


5% 
10% 


9.2% 
10.5% 


170 


5% 
10% 


8.0% 
9.0% 


200 


5% 
10% 


13.5% 
15.0% 



Table 23: Expected relative uncertainty on the measurement of quh via the process e^e~ -^ 
tiH for the combination of the 4 channels studied, for various Higgs boson masses and for 
two values of the relative uncertainty on the residual background normalisation. 
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Figure 5: The H ^ bb semileptonic channel: preselection variables (I). The signal (solid 
line) and the background (dashed line) are normalised to 1. The signal is shown for 
Mh = 120 GeV/c^. 
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Figure 6: The H ^ bb semileptonic channel: preselection variables (II). The signal (solid 
line) and the background (dashed line) are normalised to 1. The signal is shown for 
Mh = 120 GeV/c^. 
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Figure 7: The H ^ bb semileptonic channel: variables (after preselection) used for the 
neural network analysis (I). The signal (solid line) and the background (dashed line) are 
normalised to 1. The signal is shown for Mh = 120 GeV/(P. 
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Figure 8: The H ^ bb semileptonic channel: variables (after preselection) used for the 
neural network analysis (II). The signal (solid line) and the background (dashed line) are 
normalised to 1. The signal is shown for Mh = 120 GeV/cP. 
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Figure 9: The H ^ bb semileptonic channel: neural network output. Top, the signal (solid 
line) and the background (dashed line) are normalised to 1. Bottom, the signal (solid line) 
and the background (dashed line) are normalised to the expected number of events. The 
signal is shown for Mh = 120 GeV/c^ . 
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Figure 10: The H ^ bb hadronic channel: preselection variables. The signal (solid line) 
and the background (dashed line) are normalised to 1. The signal is shown for Mh = 120 
GeV/c^. 
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Figure 11: The H ^ bb hadronic channel: variables (after preselection) used for the 
neural network analysis (I). The signal (solid line) and the background (dashed line) are 
normalised to 1. The signal is shown for Mh = 120 GeV/(P. 
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Figure 12: The H ^ bb hadronic channel: variables (after preselection) used for the 
neural network analysis (II). The signal (solid line) and the background (dashed line) are 
normalised to 1. The signal is shown for Mh = 120 GeV/(? . 
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Figure 13: The H ^ bb hadronic channel: neural network output. Top, the signal (solid 
line) and the background (dashed line) are normalised to 1. Bottom, the signal (solid line) 
and the background (dashed line) are normalised to the expected number of events. The 
signal is shown for Mh = 120 GeV/c^ . 
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Figure 14: The two like sign lepton plus 6 jet channel: selection variables (I). The signal 
(solid line) and the background (dashed line) are normalised to 1. The signal is shown for 
Mh = 150 GeV/c^. 
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Figure 15: The two like sign lepton plus 6 jet channel: selection variables (II). The signal 
(solid line) and the background (dashed line) are normalised to 1. The signal is shown for 
Mh = 150 GeV/c^. 
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Figure 16: The single lepton plus 8 jet channel: preselection variables (I). The signal 
(solid line) and the background (dashed line) are normalised to 1. The signal is shown for 
Mj:/ = 150 GeV/c^. 
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Figure 17: The single lepton plus 8 jet channel: preselection variables (II). The signal 
(solid line) and the background (dashed line) are normalised to 1. The signal is shown for 
Mh = 150 GeV/c^. 
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Figure 18: The single lepton plus 8 jet channel: variables (after preselection) used for the 
neural network analysis (I). The signal (solid line) and the background (dashed line) are 
normalised to 1. The signal is shown for Mh = 150 GeV/(P. 
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Figure 19: The single lepton plus 8 jet channel: variables (after preselection) used for the 
neural network analysis (II). The signal (solid line) and the background (dashed line) are 
normalised to 1. The signal is shown for Mh = 150 GeV/(? . 



52 




« 100 

c 

a> 
o 
S 80 

E 



60 



40 



20 



1 .2 1 .4 

NN output 



- 






- 


\^^ 


1- 


: 


lM\ 


iM<i 


— 


ji"' ■ ' r 
■i.ii « 






* :\' ■ 


' "J S 'ir . ,■; 




^1 %" « 


H> ' "I'.Vl" 




^1 1- ■ 


' i ' ■ ■> 




^, .- 


' '''''i"'' 




J^ 




- 


J.lf 


i i ^'J 1' 


— 


'^" 


} # M |S 


- 




;; ;ji ;:;■ - 




m'' 


'-..'Vi'm' \ 


- 


J" 






■ ' 


*r-L - 









-0.4 



-0.2 



0.2 



0.4 



0.6 



0.8 



1.2 1.4 

NN output 



Figure 20: The single lepton plus 8 jet channel: neural network output. Top, the signal 
(solid line) and the background (dashed line) are normalised to 1. Bottom, the signal (solid 
line) and the background (dashed line) are normalised to the expected number of events. 
The signal is shown for Mh = 150 GeV/c^ . 
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Figure 21: The H ^ bb semileptonic channel: variables used in the analysis for three 
different top mass values: 170, 175 and 180 GeV/c^ (Mh = 120 GeV/c^) (I). 
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Figure 22: The H ^ bb semileptonic channel: variables used in the analysis for three 
different top mass values: 170, 175 and 180 GeV/c^ (Mh = 120 GeV/c^) (II). 
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